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Abstract 

O 

The analysis of the sum rules (for two different choices of 5-quark currents) indicate 
\q ■ that to explain experimental results it is necessary to have some one-particle lower state 

CN ■ (and pentaquark is a reliable candidate to this state) while the two-particle coupled lower 

state (NK system) should be excluded. Obtained sum rules have good stability and 
correspond to pentaquark mass mg=1.54 GeV at the appropriate choice of continuum 
threshold. Also the features of the sum rules for 2-point correlator with 5-quark currents 
are studied. This features should be taken into account to obtain sum rules correctly. 
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1 .Introduction 



Recently, the exotic narrow ( with the width less than 9MeV) baryon resonance with the 
quark content uudds and mass 1.54 GeV pQ, [2j had been discovered. It was supposed that it 
can be identified with G + baryon, predicted in 1997 by D.Diakonov, V.Petrov and M.Polyakov 
[3] in the Chiral Soliton Model as a member of antidecouplet with the hypercharge Y — 2. 

Later, the existence of this resonance was confirmed by many other groups (see for the 
review [3]). The phase analysis leads to the conclusion that the width should be even smaller 
(less than lMeV, see On the other side, in a number of experiments searches for 

this state were unsuccessful, (see jl] for the review). So one can say, that the question does 
pentaquark exist or no has yet no experimental answer. 

The explanation of small width of pentaquark was suggested in papers jHj, [HJ. It was shown, 
that if the pentaquark size is not larger than usual hadron size and its field is represented by the 
local quantum field opperator, then its decay width should be strongly suppressed by chirality 
conservation for any choice of the interpolating current. Later the numerical estimations of 
pentaquark width were also obtained in QCD sum rule approach [TQ] , [TT] , [T^J . In this paper we 
will analyze, if is it possible to predict something about the pentaquark existence using QCD 
sum rule. 

QCD sum rules (S.R.) is one of the powerful methods of analyzing such problem, and it is 
well-known, that using the sum rules for two-point correlators it became possible to predict the 
masses of usual baryoOns (see [13J). So, reasonably the question appears, is it possible to use 
the same approach for pentaquark case. 

A number of papers are devoted to this question (see, for example |E]-[IE])- Unfortunately, 
except [IB] > the analysis of the 2-point correlator in the main part of this paper was restricted at 
comparatively low dimension operators of operator product expansion (OPE) series (dimension 
d < 8). As was noted in papers [E],[IE], for the case of 5-quark current the contributions of high 
dimension operators are significant and should be taken into account. This is a consequence of 
high dimension of 5-quark current. We will discuss this in sect. 2. 

On the other point of view, such significant role of high operator contribution leads us to a 
conclusion, that the accuracy of S.R for multiquark system can not be good (particularly S.R. 
for pentaquark mass). Moreover, there is also one source of uncertainty, which is also absent in 
the case of usual hadrons. The problem is, that unlike the usual hadron case, where the masses 
of first resonance are considerably smaller than the threshold of multiparticle states( which are 
included into continuum) in the case of pentaquark the threshold of 2-particle intermediate 
state {NK, i.e nucleon and kaon) generation is even a bit smaller than a pentaquark mass. In 
this paper we will try to analyze all these problems for the pentaquark 2-point sum rules. 

The following problems will be discussed in this paper: 

1. We will discuss the properties of operator expansion for 5-quark current, and give some 
arguments, that the main contribution to S.R. give high dimension operators, but always there 
should be a strong numerical cancellation of neighboring dimensions contributions for any 
choice of current. It will be shown (and we will discuss the examples) that contribution of 
the dimensions 6 and 8 contribution cancel each other up to 90%, and the same cancellation, 
though not so strong, but also significant appears for higher dimension. 

2. It will be supposed that NK system contribution in the spectral density at the region 
close to the pentaquark mass can be interpolated by some local current, and it will be checked if 
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it is possible to saturate the 2-point sum rules only by this system plus continuum. If yes, then 
one should say, that 2-point S.R. results for 5-quark current can not be treated as an argument 
in the favor of pentaquark existence, but can be easily explained as NK system. If no, then we 
should say, that S.R. required some one-particle state in the region of masses about 1 — 2GeV 
(though the accuracy of determination of mass in any case will be not good). All the analysis 
will be done on an example of two different currents with pentaquark quantum numbers, but 
we believe that the obtained conclusion are reliable for any currents. 

2. General properties of 5-quark current correlators 

In this section we will discuss some feature of sum rules for 5-quark current correlator. 

n(p 2 ) = i J e ipx (0 | T Ve (x) Ve (0) | 0) = p ■ Tl^p 2 ) + Tl 2 (p 2 ) (1) 

Here rjg is current with quantum numbers of pentaquark and with quark content uudds 
Let us define 

<0|r7©|e> = Ai; e (2) 

where Vq is the 9 + wave function and A is the corresponding coupling constant. Then, 
using narrow resonance approximation and quark-hadron duality, one can easily write down for 
physical representation of invariant amplitudes II^p 2 and n(p 2 ) in (1) 

ImUf ys = n\ 2 5(p 2 - m 2 ) + 9{p 2 - s )ImU? CD (3) 

ImU p 2 h ys = 7i\ 2 5(p 2 - m 2 ) ■ m e + 9{p 2 - s )ImU^ CD (4) 

where ]jQ CD means the QCD calculation of the corresponding invariant amplitude, s is the 
continuum threshold and m$ is the 9 + mass. Then, equating physical representation and result 
of QCD calculation one can easily get (after Borel transformation) two sum rules for TLi (p 2 ) 
and II(p 2 ) (SRI and SR2 correspondingly). 

- I due- u / M2 ImU? CD (u) = \ 2 e- m V M2 (5) 

7T J 

- [ due- u / M2 ImIl^ CD (u) = \ 2 e~<l M2 . m e (6) 

7T J 

Here ImU® CD mean the sum of OPE series, i.e. 

ImUQ CD = J2lmU^ CD 

where ImH® D means the corresponding dimension d contribution in the operator expansion, 
(see Fig.l, where examples of diagrams for different dimension operators are given). We will 
account terms up to dimension d — 14 for SRI (eq.5) and d — 13 for SR2 (eq.6). In the 
calculation we will use the quark propagator expansion in OPE up to d = 7 operators (the first 
term of the expansion on s-quark mass will be also accounted) 
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where 



(0 | Tq a (x)q b (0) | 0) = S a \ X/ 

Z7T" L X" 

m q a5 A m*mka5 



rp . syiCX ryft o/5 ^ 

S ab -(1 - AS) - BS • 6 ab - 



(7) 



AS = + !^!o^ x 6 (g) 

96 96-24 w 

1 9 .a5/ m^x 2 6 • x 4 \ 

BS = -m s x 2 -i— 1 + —2— + (9) 

2 24 V 16 4-96-3/ V ; 

6$ = e^Gl^ejl (10) 

and 

a = -(27r) 2 (0|gg|0) 
6=(0|/G 2 |0) 

0<O|^(A B /2)G>|O) = m 2 (0|gg|0) (11) 

Here m q = m s , 5 = (0|ss|0)/(0|gg|0) for s-quark case. (It is clear, that for u, d quark 
propagators one should put m q = , and unity instead 5 ) 

We should note, that for quark propagator in the external gluon field (the last term in 
(7)) we neglect the corrections, proportional to quark mass (i.e. the term, proportional to 
m s ), because the contribution of the terms, connected with this corrections in the any order of 
OPE are smaller, than total uncertainty (due to factorization) of the main terms of the same 
dimension. Basing on the same reason we neglect all the terms, proportional to m s b and also 
some other term of high dimensions, the contribution of which was found to be smaller than 
15% of the main terms of the same dimension. (Some of such omitted terms will be discussed 
below). We neglect also all the terms, proportional to a s . That's why in the quark wave 
function expansion (we restrict ourselves to the first five terms of expansion) 

Mx) = MO) + x a V> + ^V^VV + „ V"V /3 V 7 t/;+ 

2 6 

,-y»d /-y»'"y syT 

xx x x VQy/3v 
24 

we don't account the fourth term (with three derivatives), because it contribution to S.R. 
is proportional to a s due to equation of motion, see [2D]- All terms in quark propagator, 
corresponding to quark wave function expansion (11), are accounted in BS and AS. In the 
simple case m s = the terms in the brackets in BS correspond to first, third and fifth terms in 
(11). 

Prior to we write the sum rules, let us discuss some common properties of multiquark sum 
rules. For a moment, let us for simplicity neglect the strange quark mass. Then it is clear 
from eq.(7), that for contribution of the quark condensate of the n — th power, which have the 
dimension d = 3n, there are the contributions of dimension d + 2 and d + 4 (due to second 
and third term in BS) which have fixed additional factor in comparison with the dimension d 
contribution. For example, if the contribution of d = 6 (only the part proportional to a 2 ) is 
denoted as A e 
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A 6 = x a C/x w (13) 

where C is some numerical coefficient, then d = 8 and d = 10 contribution from BS can be 
written as 

A 8 = {2m 2 /16)x a C/x 8 (14) 



rrin 2b 



no 



i * u \^ar<lJo 



•256 ' V9*' C I* (15 > 

for any current. After Fourier and Borel transformation and cancellation of continuum 
contribution with the corresponding QCD part one can find the following relations for these 
contributions 

A 8 = -2ml(E x lE 2 )As (16) 

and 

A 10 = (^+ b -)(E /E 2 )A 6 (17) 

where E n = (M 2 ) n E n , M 2 is the Borel mass square and function E n , appeared due to 
continuum extraction, is: 
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so/M 

E ^ ) = h I dzzUe ~ z (18) 



So one can see, that strong cancellation should take place for the dimension 6 and 8 con- 
tributions (A e and A 8 ), because the factor 2ml(E 1 /E 2 ) is close to 1 (at any reasonable choice 
s = 3 — 5GeV 2 and M 2 = 2 — 3GeV 2 ). Dimension 10 contribution also plays some role in this 
cancellation, (though not so significant, because A w is comparatively smaller) and should be 
also accounted. (Such situation take place for any dimension less than 10). This type cancel- 
lations really appear only for multiquark currents, for the case of usual hadrons one can easily 
see (just in the same way), that A 6 term is much more large than A$ and A w . 

Of course, one should note, that Ag and Aiq aren't the only contribution of the appropriate 
dimension. For example for d — 8 power corrections one should take into account the terms, 
when a vacuum gluon line from one quark propagator is entering in quark condensate of another 
quark line (see Fig. lb ), but all such terms are significantly smaller than As, we discuss before, 
almost for all currents (except for special cases, when A 6 = 0) and don't seriously influence on 
the cancellation. 1 

Just for the same reason for dimensions more than 10 we neglect the contribution of the 
diagrams like Fig. lb, (where two vacuum gluon lines are entering in two quark condensates) 
because their contribution is found to be smaller, than uncertainty (due to factorization hy- 
pothesis) of main terms of he same dimension (see diagrams in fig.l). It is clear, that all 
this conclusions are also valid for the terms, proportional to strange quark mass. It is easy 
to see, that this cancellation becomes negligible at large dimensions d > 11, 12, and at these 
dimensions one can safely neglect the second and third terms in BS in brackets in comparison 



1 It should be noted, that for this term there also appears a cancellation of corresponding d — 10 terms, 
which can be proved just in the same way as we have done in our example. 
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with first one. That's why the best way of calculation of the 2-point correlator for 5-quark 
current is to extend OPE at least to dimension 12,13 to avoid the uncertainty due to possible 
cancellations with higher dimensions. 

Another feature of multiquark currents is domination of the high dimension correlators (see 
jH],PHl), even in spite the cancellation we discuss before. The reason is that lower dimension 
operators are strongly suppressed because of continuum cancellation (and also due to small 
numerical factor appearing after Fourier transformation), which lead to the following factor 
(for SRI, for example) 

95-2i 

(z-1)! V ; 

Here i = 8—d/2, and a small factor in front of E appears due to Fourier transformation. One 
can see that for low dimension operators (d < 6) the function E n , appearing due to continuum 
cancellation, is extremely small (it totally compensate the large value of M 2n ). This reflects the 
fact, that for the multiquark system the perturbative contribution (as well as low dimension 
operators) mainly corresponds to the multiparticle states (the region of quark-hadron duality), 
and not to the resonance region. 2 

Thus we can conclude that for multiquark current it is important (unlike the usual hadron 
case) : 

1. to take into account high dimensions up to dimension 12,13, 

2. if the OPE series is terminated by dimension less than 12, keep only those terms for 
which the cancellation terms of the next order are accounted - in other case the answer can be 
changed drastically. 

3. The pentaquark mass is determined from the ratio of eq.'s (6) and (5). As was shown in 
Ref. i n the case of multiquark baryon currents the OPE series convergence is much better 
in this ratio, than in (5) or (6). (The QCD sum rules for baryons at large N c were considered in 
[19J.) So, one may believe, that the values of the pentaquark masses found below are reliable, 
if the formulated above hypothesis about the representation of pentaquark field by local field 
operator is valid. 



3. Sum rules 



Let us consider now the sum rules (both chirality conserving and violating, - SRI and SR2 
correspondingly, eq.'s (2,3)) for two following currents with quantum number of 9 + pentaquark 

J A = e abc e def e 9cf (u aT Cd b )(u dT C^ 5 d e )^cs g (19) 

J B = {e abc [(d a C^d b )u c - (d a C^u b )d c ] ■ s 7/ ,7 5 m + u <-> d}/V2 (20) 

The reason of such choice of currents as characteristic examples for 9 + sum rules is that 
the structure of this currents corresponds to two most reliable cases: the first - to 3-quark 
system( analogous to baryon) multiplied to .fT-meson-like part and the second - two diquarks 

2 We would like to note, that even at a special current when contribution of d=6 due to quark condensate 
square is zero, the contributions of high dimensions are found to be large enough. 
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and antiquark. The large number of similar currents(though not just the same) were obtained 
in a number of papers (see, for example ^3]- [TS]). As we have discussed above, the role of the 
high dimension corrections is very significant. So we account power corrections up to dimension 
14,13 (with restrictions we discuss in the previous section). Then the final sum rules SR1,SR2 
can be written in the following form (we will further mark them SR1A and SR2A for the current 
j a and SR2B and SR2B for the current 



M 2 K{ A > B) E 5 + R^E 3 + R[ A > B) E 2 + R^E, + R[ A ' B) E + 



+ R {a ' b) /M 2 + r£ B) /M 4 = AV 



rrv 



J/M 2 



(21) 



M 2 K (A,B) R (A,B)z 5 + R (AB)^ + R ^,B)^ + ^AB)^ + ^B)^ 



+ R[f B) E + r[ A > B) /M 2 = m e \ 2 e-<l u 



(22) 



where 







1925 




24 





2625 



192 







-2aS 




Rf = 40m s a 2 + —abd 
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Rg = 40a 3 8 - 50m s mla 2 - ^«.« () 

8 



m n ab5 



R 



325 



ii 



m a 5 + 10m m s a 



pa 
^13 



185 
1}6~ 



12 

k? = ^ 

1 35 



160 



-m,a 



18 



i?f = !(^-|m.a(l + 2*)) 



35 

is' 



R% = 2l(3 m l msa s _ 3a 2(! + 5 )) 
18 

^f = |« 2 ^(14 + 3o) 
33 



( -6m s a 3 (2 - 5) - —m A Q a 2 + a 2 o(2 + 5/2) 
| (-8(1 + o)a 4 + m^a 3 (ll - ^ - m*a 2 &(^ + ^ 



m 2 a 4 (^ + y + mjra s fl 3 (-1 + ^ + mfc 2 6(3/8 + 3/165)) 



i2f 



1? 



B = _ {1 _ 6)m 2 a 



R 



13 



— (12(1 - 5/2)m s a + (— + 35)b)a 
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27 



—69 9 9 15 

12(1 + 5)a 2 + (— + -5)m s m 2 a ~(^ + ^)m 2 b ) a 



R?i = ^(- 2 ( 6 + ±75)m 2 a 2 + (^ + ^m.m^a 

5 / 3 2"i \ 

& = 27 (-4(8 + o)m s a 3 + (- + y o)m 4 a 2 + 2a 2 oJa 



Here the indexes A,B denote the sum rules for the currents j a , jb (see eq. (19,20)), R n 
means the contribution of the operators of dimension n, and A 2 = (47r) 8 A 2 , where A is defined 
in (2). 

Continuum contributions are accounted in the left-hand sides (l.h.s) of the s.r. resulting in 
appearance of the factors E n . 

We can easily see that the strong cancellation we discuss in the previous section, really 
appears (for example R 6 and R$ e.t.c). 

The values of A 2 , determined from eqs.(21) and (22) with the continuum threshold chosen 
as So = 3.5GeV 2 are plotted in Fig. 2a for the value of mg = 1.54GeV, and the value of mg 
obtained as a ratio of (22) to (21) is plotted in Fig.2b. The parameters were taken in accord with 
the recent determination of QCD condensates [21], [22] at normalization point /i 2 = 2GeV 2 : 
a = 0.63Ge^ 3 , b = 0.2AGeV\ m 2 = 0.8GeV 2 , m s = Q.lbGeV, 5 = 0.8. 

We vary the value of the continuum threshold within the range 3-4.5GeV 2 . One can easily 
see, that sum rules dependence on the Borel mass are rather stable (both SRI and SR2) and 
also the value of the pentaquark mass is also stable enough. On the other side, the result 
strongly depends on the continuum threshold value. This fact is not unexpected, such strong 
dependence on the continuum threshold reflects the fact we have discussed before: in the 
pentaquark case the physical threshold of the first multiparticle state (nucleon plus kaon) in 
physical representation in physical representation of 2-point correlator, (which usually is much 
larger then first resonance and is included in the continuum), is even a little smaller than the 
mass of the first resonance (pentaquark 6 + ) mass. That's why in the pentaquark case the 
"continuum threshold" is just the free parameter, which can not be directly associated with 
higher resonance masses or physical multiparticle states threshold, and it will be fixed if one 
demands that the ratio of the 2-point sum rules SR2/SR1 should be equal to the experimental 
value. That's why, it seems, that it is principally impossible to predict the pentaquark mass 
(unlike the case of usual hadrons)from 2-point sum rule for the pentaquark case. The only thing 
we can find from 2-point sum rule is: using the experimentally known value of mass, choose 
appropriate currents, which give reasonable sum rules (SRI and SR2) and fix the such value 
of the parameter sq that the SR2/SR1 ratio becomes about pentaquark mass -if it is possible 
to do. (Of course, as soon as this parameter is fixed, it will be used in any other sum rule 
with this current, according to the usual logic of the sum rule approach). So we see, that for 
our case we can conclude that both currents j a jb can be treated as candidates to pentaquark 
current, because they give sum rule with a good stability and agree with experimental value of 
the pentaquark mass at so about 3.5GeV 2 , but one can't say, that sum rules (21,22) themselves 
predict the pentaquark existence. 

And yet, surprisingly, it was found to be possible to get more information from 2-point 
correlator and get some predictions about pentaquark. We will discuss this in the next section. 

4. The estimation to the KN contribution to the sum rules 

The main question we try to answer in this section is: is it possible to explain the exper- 
imental results as some coupled baryon-meson system without any one-particle intermediate 
state. Let us suppose now, that the size of the NK system (nucleon and kaon), is not too large 
(less than 0.5fm), so it can be extrapolated by some interpolating current (just the same as for 
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pentaquark). To account the fact that this NK system is not local let us define 



(0 | 77 I nk) = (aP + /3m)7 5 tv, 



(25) 



where a, f3 are two unknown formfactors, P is the four- moment, carried by nuclon-kaon 
system, and v n is nucleon spinor and m is nucleon mass. It is clear, that if we will be interested 
in the region of the order of the pentaquark mass 1.54GeV, which is close to (m n + m k ) 2 , (where 
m n and m k are nucleon and kaon masses correspondingly), one can suppose that formfactors 
are more or less constant in this region. 

Then, substituting (25) into the relation ImU. — X)(0 | 77 | nk){nk \ f) | 0) + cont 
and integrating of the all intermediate states momenta, one can easily obtain for physical 
part of sum rules 



where 



ImU = Ppi(u) + mp 2 (u) + cont.contr 



u 



(1 — z) 2 \ ml z 2 m 2 (m 2 —ml) 



+ 



2u 



(26) 



(27) 



p 2 = —en ((z — l)m(u — zm ) — zm ■ m k ) 

07T 



(28) 



and z = (3/a, c = \J (1 — {jri^ + m\)/u) 2 — (2m n mk/u) 2 

Then after Borel transformation one can easily write sum rules for both kinematical struc- 
tures in the form, analogous to (21,22), where all notations of the left QCD side were defined 
in (23,24). 



(A,B) 



E 5 + R[ A ^E 3 + R[ A ' B) E 2 + R^E, + R[ A ' B) E Q + 



+ R[ A ' B) /M 2 + R[f/M 



(A,B) 



M 2 
(4^f 



so 

= J due- u ' M2 Pl /Ti 



(mn+mk) 2 



K. 



(A,B) 



T?(A,B)p R(AB)p p(AB)p p(AB)-p o(AB)p 



+ r[ a ' b) e + r[ a ' b) /m 2 



M 2 



so 

= J due- u ' M2 p 2 /n 



(29) 



(30) 



(m n +m k ) 2 

On the other side at the pure SU3f chiral limit ( when the mass of kaon is zero) and the 
momentum of kaon q = the matrix element (25) should tend to zero. It can be easily shown 
using the equation of motion, that then for formfactors we have the relation a — f3. Because 
in our case we are interesting at physical region about pentaquark mass, which is close to NK 
system generation threshold [m n + m k ) 2 , then all momentum of particles are rather small and 
it is reasonable to demand that formfactors can differ at the same order as SU3f violation, i.e. 
about 30%. 
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So for each current we obtain two different sum rules, and the question is: is it possible to 
find such values of formfactors (more or less of the same value, i.e z vary at 0.7 to 1.3, and a 
is free), which satisfy both sum rules (29,30). (One should note, that the parameter so in this 
case is not the same as for the sum rules for the pentaquark case; we will vary it, but, of course, 
not very far from the pentaquark mass square region). The numerical analysis shows, that for 
both currents (19,20) it found to be impossible to satisfy sum rules (29) and (30) at the same 
time for all reliable values of z. For example in Fig.3a,b, the value £, 

£ = l/(4Va) 2 (31) 

obtained from SRI (29) and SR2 (30)are shown for both choices of current j A (Fig. 3a), 
and j B (Fig.3b), for the case z = 1.3. One can see, that results of SRI and SR2 totally 
contradict each other. For the smaller values of z the contradiction becomes even more sronger 
(for example at z = 1 the results of sum rules (29,30) strongly differ by modulo and moreover, 
have the opposite sign). So we can conclude, that the sum rules for the currents we consider 
indicate that to explain experimental results it is necessary to have some one-particle lower 
state (and pentaquark is a reliable candidate to this state) while the two-particle coupled lower 
state (NK system) should be excluded . 

4. Summary 

The main results of this paper are the following: 

1. The features of the sum rules for 2-point correlator with 5-quark currents are studied. It is 
shown, that on one side, the contributions of high dimension operators are significant, and, on 
the other side, the strong cancellation of the neighboring dimensions contribution should take 
place for any choice of the interpolating current. These features should be taken into account 
to obtain sum rules correctly, as is discussed in sect. 2. The best way, from our point of view, is 
to take into account OPE series up to dimension 13,14, where one can neglect the cancellations. 

2. Analysis of the sum rules for 2-point correlator with 5-quark currents for two different 
currents show, that it is impossible to explain the experimental result on pentaquark, supposing 
the existence only 2-particle lower state (NK -system) in the physical representation of sum 
rules. For this reason one should suppose the existence of any one-particle lower state (like 
pentaquark) in this mass region. Obtained sum rules (21,22) have good stability and correspond 
to pentaquark mass me— 1-54 GeV at the appropriate choice of continuum threshold. 

It is very necessary to note, that all conclusions and results are obtained without a possible 
instanton contribution. On the other side, it is well known, that the instanton contribution can 
be significant, especially at large dimension of OPE. (see [22] - |25j ) 

That's why the instanton contribution can influence the obtained results and conclusions 
and we are planning to investigate this problem in future. 

Author thanks B.L. Ioffe for many significant advises and useful discussions. 
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Figure 1: Examples of the diagrams, corresponding to the different operators contribution. 
Circles mean derivatives 
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Figure 2: a)Sum rules (eqs.21,22) for currents j A (1A and 2A correspondingly) and j b (1B,2B). 
b) the pentaquark mass obtained as ratio eq.22/eq.21 (marked SRA)for currents j A and the 
same for currents j B 
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Figure 3: a) The value of £ (see eq.31) obtained from sum rules eq.29 (SR1A,) and eq.30 (SR2A) 
for current j A . b) the same for currents j B 
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